Continuous variable entanglement distillation of Non-Gaussian Mixed States 
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Many different quantum information communication protocols such as teleportation, dense cod- 
ing and entanglement based quantum key distribution are based on the faithful transmission of 
entanglement between distant location in an optical network. The distribution of entanglement in 
such a network is however hampered by loss and noise that is inherent in all practical quantum 
channels. Thus, to enable faithful transmission one must resort to the protocol of entanglement 
distillation. In this paper we present a detailed theoretical analysis and an experimental realization 
of continuous variable entanglement distillation in a channel that is inflicted by different kinds of 
non-Gaussian noise. The continuous variable entangled states are generated by exploiting the third 
order non-linearity in optical fibers, and the states are sent through a free-space laboratory channel 
in which the losses are altered to simulate a free-space atmospheric channel with varying losses. 
We use linear optical components, homodyne measurements and classical communication to distill 
the entanglement, and we find that by using this method the entanglement can be probabilistically 
increased for some specific no n- Gaussian noise channels. 



PACS numbers: 42.50.Lc,42.50.Dv,42.81.Dp,42.65.Dr 

I. INTRODUCTION 

Quantum communication is a promising platform for 
sending secret messages with absolute security and de- 
veloping new low energy optical communication sys- 
tems p. Such quantum communication protocols re- 
quire the faithful transmission of pure quantum states 
over very long distances. Heretofore, significant experi- 
mental progress has been achieved in free space and fiber 
based quantum cryptography where communication over 
more than 100 km have been demonstrated [2H4]. The 
implementation of quantum communication systems over 
even larger distances - as will be the case for transatlantic 
or deep space communication - can be carried out by us- 
ing quantum teleportation. However, it requires that the 
two communicating parties share highly entangled states. 
One is therefore faced with the technologically difficult 
problem of distributing highly entangled states over long 
distances. The most serious problem in such a transmis- 
sion is the unavoidable coupling with the environment 
which leads to losses and decoherence of the entangled 
states. 

Losses and decoherence can be overcome by the use 
of entanglement distillation, which is the protocol of ex- 
tracting from an ensemble of less entangled states a sub- 
set of states with a higher degree of entanglement [6]. 



Distillation is therefore a purifying protocol that selects 
highly entangled pure states from a mixture that is a 
result of noisy transmission. This protocol has been ex- 
perimentally demonstrated for qubit systems exploiting 
a posteriori generated polarization entangled states [7]- 
[TT] . Common for these implementations of entanglement 
distillation is their relative experimental simplicity; only 
simple linear optical components such as beam splitters 
and phase shifters are used to recover the entanglement. 
This inherent simplicity of the distillation setups arises 
from the non- Gaussian nature of the Wigner function of 
the entangled states. It has however been proved that 
in case the Wigner functions of the entangled states are 
Gaussian, entanglement distillation cannot be performed 
by linear optical components, homodyne detection and 
classical communication [T2HT4] . This is a very impor- 
tant result since it tells us that standard continuous vari- 
able entanglement generated by e.g. a second-order or 
a small third-order non-linearity cannot be distilled by 
simple means as these kinds of entangled states are de- 
scribed by Gaussian Wigner functions. 

Several avenues around the no-go distillation theo- 
rem have been proposed. The first idea to increase the 
amount of CV entanglement was put forward by Opatrny 
et al. [15 . They suggested to subtract a single photon 
from each of the modes of a two-mode squeezed state 
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using weakly reflecting beam splitters and single pho- 
ton counters, and thereby conditionally prepare a non- 
Gaussian state which eventually could increase the fi- 
delity of CV quantum teleportation. This protocol was 
further elaborated upon by Cochrane et al. [16 and Oli- 
vares et al. [T?] . Other approaches relying on strong cross 
Kerr nonlinearities were suggested by Duan et al. [19] [20] 
and Fiurasek et al. [21]. The usage of such non-Gaussian 
operations results in non- Gaussian entangled states. To 
get back to the Gaussian regime, it has been suggested to 
use a Gaussification protocol based on simple linear op- 
tical components and vacuum projective measurements 
(which can be implemented by either avalanche photodi- 
odes or homodyne detection) [22] . Distillation including 
the Gaussification protocol was first considered for pure 
states by Browne et al. [22] but later extended to the 
more relevant case of mixed states by Eisert et al. [23 . 

Due to the experimental complexity of the above 
mentioned proposals, the experimental demonstration of 
Gaussian entanglement distillation has remained a chal- 
lenge. A first step towards the demonstration of Gaussian 
entanglement distillation was presented in ref. [24] where 
a modified version of the scheme by Opatrny et al. [T5] 
was implemented: single photons were subtracted from 
one of the two modes of a Gaussian entangled state us- 
ing a nonlocal joint measurement and as a result, an in- 
crease of entanglement was observed. Recently, the full 
scheme of Opart ny et al. was demonstrated by Taka- 
hashi et al. [25]. They observed a gain of entanglement 
by means of conditional local subtraction of a single pho- 
ton or two photons from a two-mode Gaussian state. 
Furthermore they confirmed that two-photon subtraction 
also improves Gaussian-like entanglement. 

In the work mentioned above, only Gaussian noise has 
been considered as for example associated with constant 
attenuation. Gaussian noise is however not the only kind 
of noise occurring in information channels. E.g. if the 
magnitude or phase of the transmission coefficient of a 
channel is fluctuating, the resulting transmitted state 
is a non- Gaussian mixed state. Because of the non- 
Gaussianity of the transmitted state, the aforementioned 
no-go theorem does not apply and thus Gaussian trans- 
formations suffice to distill the state. Actually, such a 
non-Gaussian mixture of Gaussian states can be distilled 
and Gaussified using an approach [26] related to the one 
suggested in ref. [22]. Alternatively, it is also possible to 
distill and Gaussify non- Gaussian states using a simpler 
single-copy scheme which is not relying on interference 
but is based on a weak measurement of the corrupted 
states and heralding of the remaining state [27]. Such 
an approach has been also suggested for cat state purifi- 
cation [28], coherent state purification [29] and squeezed 
state distillation [30] . 

The distillation of Gaussian entanglement corrupted 
by non- Gaussian noise was recently experimentally 
demonstrated in two different laboratories. More specifi- 
cally, it was demonstrated that by employing simple lin- 
ear optical components, homodyne detection and feedfor- 



ward, it is possible to extract more entanglement out of a 
less entangled state that has been affected by attenuation 
noise [27 or phase noise [3T] , 

In this paper we elaborate on the work in ref. [27] , 
largely extending the theoretical discussion on the char- 
acterization of non-Gaussian entanglement and, on the 
experimental side, testing our distillation protocol in new 
attenuation channels. 

The paper is organized as follows: In Section [TTJ the 
entanglement distillation protocol utilized in our experi- 
ment is fully discussed. In Section [TTT| the experimental 
setup for realization of the entanglement distillation is 
described, and the experimental results are shown and 
discussed in Section HVl 



II. THE ENTANGLEMENT DISTILLATION 
PROTOCOL 

The basic scheme of entanglement distillation is il- 
lustrated in Fig. [2] The primary goal is to efficiently 
distribute bipartite entanglement between two sites in 
a communication network to be used for e.g. teleporta- 
tion or quantum key distribution. Suppose the two-mode 
entangled state (also known as Einstein-Podolsky-Rosen 
state (EPR)) is produced at site A. One of the modes is 
kept at A's site while the second mode is sent through 
a noisy quantum channel. As a result of this noise, the 
entangled state will be corrupted and the entanglement 
is degraded. The idea is then to recover the entangle- 
ment using local operations at the two sites and classical 
communication between the sites. To enable distillation, 
it is however required to generate and subsequently dis- 
tribute a large ensemble of highly entangled states. After 
transmission, the ensemble transforms into a set of less 
entangled states from which one can distill out a smaller 
set of higher entangled states. 

A notable difference between our distillation approach 
and the schemes proposed in Refs. [22] [26] is that our 
procedure relies on single copies of distributed entangled 
states whereas the protocols in Ref. [22j [26] are based 
on at least two copies. The multi copy approach relies 
on very precise interference between the copies, thus ren- 
dering this protocol rather difficult. One disadvantage 
of the single copy approach is the fact that the entan- 
gled state is inevitably polluted with a small amount of 
vacuum noise in the distillation machine. This pollution 
can, however, be reduced if one is willing to trade it for 
a lower success rate. 

Before describing the details of the experimental 
demonstration, we wish to address the question on how 
to evaluate the protocol. The entanglement after distil- 
lation must be appropriately evaluated and shown to be 
larger than the entanglement before distillation to en- 
sure a successful demonstration. One way of verifying 
the success of distillation is to fully characterise the input 
and output states using quantum tomography and then 
subsequently calculate an entanglement monotone such 
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as the logarithmic negativity. However, in the experi- 
ment presented in this paper (as well as many other ex- 
periments on continuous variable entanglement) we only 
measured the covariance matrix as such measurements 
are easier to implement. The question that we would 
like to address in the following is whether it is possible 
to verify the success of distillation based on the covari- 
ance matrix of a non-Gaussian state. 



A. Entanglement evaluation 

In order to quantify the performance of the distillation 
protocol, the amount of distillable entanglement before 
and after distillation ought to be computed. It is how- 
ever not known how to quantify the degree of distillable 
entanglement of non-Gaussian mixed states [6j [32j [33] . 
Therefore, as an alternative to the quantification of the 
distillation protocol, one could try to estimate qualita- 
tively whether distillation has taken place by compar- 
ing computable bounds on distillable entanglement be- 
fore and after distillation. First we will have a closer 
look at such bounds. 



1. Upper and lower bounds on distillable entanglement 

Although it is unknown how to find the amount of dis- 
tillable entanglement of non-Gaussian mixed states, we 
can easily find the upper and lower bounds by comput- 
ing the logarithmic negativity and the conditional en- 
tropy, respectively [34-36 . These bounds can be found 
before and after distillation, and the success of the dis- 
tillation protocol can be unambiguously proved by com- 
paring these entanglement intervals: If the entanglement 
interval is shifted towards higher entanglement and is not 
overlapping with the interval before the distillation, the 
distillation has proved successful. In other words, distil- 
lation has been performed if the lower bound after the 
protocol is larger than the upper bound before. This is 
illustrated in Fig. [TJa). 

It has been proved that for any state, the log-negativity, 



LN(p) = log 2 (2A/> 1) = log 2 \\p T % 



(1) 



is an upper bound on the distillable entanglement; Ed < 
LN(p) [34]. Here p is the density matrix of the state, 
||p TA || is the trace norm of the partial transpose of the 
state with respect to subsystem A, and the negativity is 
defined as 



Af(p) 



n T A\ 



(2) 



The negativity corresponds to the absolute value of the 
sum of the negative eigenvalues of p TA and it vanishes for 
non-entangled states. 

In our experiment we were not able to measure the 
density matrix and thus compute the exact value of the 



negativity. We therefore use another (more strict) upper 
bound that is experimentally easier to estimate. As the 
negativity is a convex function we have 



(3) 



where pi denotes the zth hermitian component in the 
mixed state, and pi is the weight for the ith component 
with pi > and YliPi = 1- Using this result we can find 
an upper bound on the log-negativity for mixed states: 



LN^PiPi) < log 2 ^ + 2j2pMpi)J ■ 



(4) 



This upper bound for the log-negativity will later be used 
to compute an upper bound for the distillable entangle- 
ment. 

Another entanglement monotone is the conditional en- 
tropy. In contrast to the log-negativity, the conditional 
entropy yields a lower bound on the distillable entangle- 
ment: Ed > S(pa) — S(p) [23 [36], where p is the density 
matrix corresponding to Gaussian approximation of the 
state and pA is the reduced density matrix with respect 
to system A. The entropies of the states can be calculated 
from the covariance matrix, CM, using 

5(p A ) = /(detA), 



1 x- 
— log 2 ( — 



*i i x - 
-log 2 ( — 



: ), (5) 



where 



Ml,2 



/ 7 ± ^ -4detCM 



(6) 



are the symplectic eigenvalues of the covariance density 
matrix and 7 = det A + detB + 2detC. Here A, B, 
and C are submatrices of the covariance matrix: CM = 
{A,C;C T ,B}. 

It is important to note that this lower bound is very 
sensitive to excess noise of the two-mode squeezed state. 
Even for a small amount of excess noise, the lower bound 
approaches zero and thus is not very useful. This is il- 
lustrated in FigjTJb) which shows the distillable entan- 
glement intervals before and after distillation of a noisy 
entangled state. Although the distillation protocol might 
remove the non- Gaussian noise of the state, the Gaussian 
noise of the state persists, and thus the entropy (that is 
the lower bound on distillable entanglement) will remain 
very low even after distillation. This results in an over- 
lap between the two entanglement intervals and thus the 
comparison of computable entanglement bounds fails to 
witness the action of distillation in terms of distillable 
entanglement. 
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(a) distillation with initially pure state 
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sified in the distillation process and the approach cannot 
be applied. For such cases, however, we will resort to 
evaluations of the Gaussian part of the state in terms of 
Gaussian entanglement. 



Entropy 
Lower bound 



Log-Neg 
Upper bound 



Distillable 
entanglement 



(b) distillation with initially mixed states 
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FIG. 1. Schematic demonstration of entanglement distillation 
of non-Gaussian mixed states. In figure (a), the distillation 
with a pure state is illustrated via the shift of the entangle- 
ment interval composed by the upper and lower bounds on 
distillable entanglement before and after the distillation pro- 
tocol. Figure (b) shows the distillation with mixed states, the 
lower bound of which does not manifest increase even for a 
small excess noise in the state. 



3. Gaussian entanglement 

In addition to an increase in distillable entanglement 
and logarithmic negativity, the protocol can also be eval- 
uated in terms of its Gaussian entanglement. Although 
the Gaussian entanglement is not accounting for the en- 
tanglement of the entire state (but only considers the 
second moments), it is quite useful as it directly yields 
the amount of entanglement useful for Gaussian proto- 
cols, a prominent example being teleportation of Gaus- 
sian states. 

In a Gaussian approximation, the state can be de- 
scribed by the covariance matrix CM [37] . The loga- 
rithmic negativity (LN) can then be found as 



LN 



■ log 2 Vrr 



(7) 



where v m in is the smallest symplectic eigenvalue of the 
partial transposed covariance matrix. The symplectic 
eigenvalues can be calculated from the covariance matrix 
using 



^1,2 



'S± VS 2 -4det CM 



(8) 



2. Logarithmic negativity 

In our experiment, the entangled states possess a large 
amount of Gaussian excess noise and thus the prescribed 
method is insufficient to prove the act of entanglement 
distillation using distillable entanglement as a measure. 
However, in certain cases we can use the logarithmic neg- 
ativity as a measure to witness the act of entanglement 
distillation even though we only have access to the co- 
variance matrix as we will explain in the following. 

First we note that in general, the Gaussian logarith- 
mic negativity is an insufficient measure of entanglement 
distillation of non-Gaussian states as this measure only 
yields an upper bound, and with upper bounds of LN 
both before and after distillation a conclusion cannot be 
drawn. However, if the state after distillation is perfectly 
Gaussified its Gaussian LN becomes the exact LN, and if 
this exact value of LN is larger than the upper bound of 
LN before distillation (computed from (eqn.|4|), one may 
successfully prove the action of entanglement distillation 
entirely from the covariances matrices. This condition 
will be used for some of the experiments presented in 
this paper. More specifically, we will use this approach 
for testing entanglement distillation in a binary trans- 
mission channel. For other transmission channels inves- 
tigated in this paper, the state will not be perfectly Gaus- 



where 5 = det A + det B — 2 det C, A, B and C represent 
the submatrices in the correlation matrix [34]. Then by 
finding the smallest eigenvalue of the covariance matrix 
and inserting it in eqn. ([7|, a measure of the Gaussian 
entanglement of the state can be found. 



B. Theory of our protocol 



Beam A 

^ Conditioning — 



Lossy channel 



Conditioning - 




Beam B 



FIG. 2. Schematics of the entanglement distillation protocol. 
A weak measurement on beam B is diagnosing the state and 
subsequently used to herald the highly entangled components 
of the state. 

We now undertake our experimental setup a theoretical 
treatment in light of the results of the previous section. 

The schematic of our protocol is shown in Fig. [2] The 
two-mode squeezed or entangled state is produced by 
mixing two squeezed Gaussian states at a beam splitter. 
The squeezed states are assumed to be identical with vari- 
ances Vs and Va along the squeezed and anti-squeezed 
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quadratures, respectively. The beam splitter has a trans- 
mittivity of T$ and a reflectivity of Rs = 1 — T5. One 
mode (beam A) from the entangled pair is given to Alice 
and the other part (beam B) is transmitted through a 
fading channel. The loss in the fading channel is char- 
acterized by the transmission factor < rj(t) < 1 which 
fluctuates randomly. The probability distribution of the 
fluctuating attenuation can be divided into N different 
slots each associated with a sub-channel with a constant 
attenuation. The transmission of sub-channel i is rji and 
it occurs with the probability pi so that Y^=iPi = 1. 
For a particular ith sub-channel with transmission of rji, 
the transmitted state is Gaussian and can be fully char- 
acterised by the covariance matrix CM^: 



CM, 
B 



CT B ? 



( v APii 
V v BP ,i 



c = 



V AX ,i o 
V BX ^ 

C x ,i 

Cpa 



.(9) 



From the Wigner function the second moments of the 
quadratures can be calculated through integration: 

(ZY^ = J dx A dx B dp A dp B zyW(x A ,x B ,p A ,p B ) 
= ^Pi / dx A dx B dp A dp B zyWi(x A ,x B ,p A ,p B )i\A) 

i ^ 

where Z, Y = X A ,X B ,P A ,P B . As the first moments 
of the vacuum squeezed states in both the quadratures 
are zero, the variances directly correspond to the second 
moments. Therefore, the elements of the total covariance 
matrix are simply the convex sum of the symmetrical 
moments (10): 



(ZY) = ^2pi(ZY)i 



(15) 



Since all the moments (10) are just linear combinations 

the covariance 



where the elements are given by: 



of the transmission factors rji and y/rjl, 
matrix of the mixed state has the following elements: 



Vaxa 


= T s Vs + R s V A , 




V AX 


= T s Vs + R s V A , 




V BX ,i 


= ViiTsVA + R S V S ) + (1 - 


m), 


Vbx 


= (v) (TsV A + RsV s ) + (1 - 


(v)), 


Vap,i 


= T s V A + R s Vs, 




V A p 


= T S V A + R S V S , 




V BP ,i 


= Vi(T S V s + R S V A ) + (1 - 


Vi), 


V bp 


= (v) (TsVs + R S V A ) + (1 - 


(v)), 


C x ,i 


= -C P ,i = VmVRsT s (v A 


-V s ). (10) 


C x 


= -C P = <V5t> VRsT s (V a - 


-V s ) 



Then according to eqn. Q we can find an upper bound 
for the log-negativity of the state after transmission in 
the fluctuating channel (using||p^|| = l/^ m i n ,i) : 



LNQ^piPi) <log 2 J>*A 



(11) 



where Vmin,i corresponds to the smallest symplectic 
eigenvalue of the zth partial transposed covariance ma- 
trix. This means that the right hand side of this ex- 
pression is also an upper bound on the distillable en- 
tanglement of the non-Gaussian noisy state. Therefore, 
to truly prove that the entanglement has increased, this 
bound must in principle be surpassed. 

We now consider the Gaussian entanglement of our 
states using the Wigner function formalism. The Wigner 
function of the total state and the zth state can be de- 
scribed as 

N 



W(X,P) = ^>^(X,P), 



Wi(X,P) 



exp 



(-xv^x^ 



PVpJpr) 



(12) 



where X = (x Ai xb\ 
are given by 



47r 2 v /detV x ,i detV P>i 
and P = (p a ,Pb)- Vx,i and \ P , 



V A x,i Cx,i 

Cx,i V B X,i 



P,i 



V A p,i 

Cp,i 



Cp y i 
Vpp,i 



13) 



(16) 

where the symbol (•) denotes averaging over the fluctu- 
ating attenuations. Comparing this set of equations with 
the set in (10) associated with the second moments for 



the single sub-channels, we see that the attenuation co- 
efficient r) is replaced by the averaged attenuation (77), 
and y/fj is replaced by (s/rj). It is interesting to note 
that if the attenuation factor is constant (which means 
that the transmitted state will remain Gaussian) there 
will always be some, although small, amount of Gaus- 
sian entanglement left in the state. On the other hand, 
if the attenuation factor is statistically fluctuating as in 
our case, the Gaussian entanglement of the non-Gaussian 
state will rapidly degrade and eventually completely dis- 
appear. 

To implement entanglement distillation, a part of the 
beam B is extracted by a tap beam splitter with trans- 
mitt ivity T. A single quadrature is measured (for ex- 
ample the amplitude quadrature, X t ) and based on the 
measurement outcome the remaining state is probabilis- 
tically heralded; it is either kept or discarded depending 
on whether the measurement outcome is above or below 
the threshold value x t h- The conditioned Wigner func- 
tion of the output signal state after the distillation is 

W p {x A ,p A ,x' B ,p B ) = 

/ dx t I dp t ^2piWi(x Al p A ,x B ,PB)Wo(x Vl p v ). 

J x t h J-00 i=1 

(17) 
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where x B = VTx' B - y/l^Tx u p B = y/Tp' B + V^Tpt , 
x v = VTx t - y/l - Tx' B and p v = y/l - Tp t + y/Tp' B , 
the Wigner function Wo(x v ,p v ) represents the vacuum 
mode entering the asymmetric tap beam splitter. After 
integration, the Wigner function can be written as 

W p (x A ,PA,x' B ,p' B ) = 

—-y^PiW'x^XA^x^xth) x W'p^paiPb)- (18) 
Ps 

This is a product mixture of two non-Gaussian states 
which should be compared to the state before distillation 
which was a mixture of Gaussian states. Ps is the total 
probability of success. 

The X related elements of the covariance matrix can 
be calculated from this Wigner function directly by com- 
puting the symmetrically ordered moments: 



with 



(Xa)? 



{X' b )\ 



(x B 2 )r = 



(X A X B )[ 



(Xa) p 
{X' b ) P 
(X 2 A ) P 

(x B Y 

(X A X' B ) P 



EMx' A )r 

Ps ' 
Ps ' 

Y,Mx'2)f 

Ps ' 
Ps ' 

Y,iPi{x A x' B )f 

Ps 



(19) 



: exp 



°th 



fTR^VBXi ~ 1) 



exp 



h th 



RC\ flu 



: exp 



°th 



V AX ,i 



Erfc 



Xth 



RT(V' - lfx th 



BT(Vi 



BXA 



■ exp 



Vbxa 



b th 



Erfc 



x t h 



VTR(V{, Xl -l)C Xl 



TCx,iV BXi 



Erfc 



exp 



xth 



b th 



^DX, 



2V, 



DX,i 



(20) 



where V BX i = RV BX ,i 
detected mode and 



T is the output variance of the 



PSA 



1 



Erfc 



Xth 



(21) 



is the success probability of distilling the l-th constituent 
of the mixed state. The total probability of success Ps 
is then given by P s = ^ZiPiPs,i- 

Since the first moments of the P quadrature are van- 
ishingly small, the P related elements of the covariances 
matrix are directly given by 

i 

(PbY = ^- TpiPs,i(TVBP,i + R), 



(PaP'b) 



i \P 



T 

~P~s~ 



(22) 



The covariance matrix CM P can then be constructed 
from these elements. This covariance matrix fully char- 
acterizes the Gaussian part of the state and thus yields 
the Gaussian log- negativity by using eqn. 

As we discussed in Section [II A[ to successfully demon- 
strate entanglement distillation of non-Gaussian states, 
the upper bound on distillable entanglement before dis- 
tillation must be surpassed by the lower bound on distil- 
lable entanglement after the distillation (see also Fig. [I]). 
Due to the fragility of the lower bound, this can be only 
achieved for almost pure states as mentioned above. A 
theoretical demonstration is given in Fig. [3] Here we 
consider the transmission of entanglement in a channel 
which is randomly blocked: The entangled state is per- 
fectly transmitted with the probability p\ = 0.2 and 
completely erased with the probability p2 = 0.8. We 
assume the two squeezed states which produce entan- 
glement have variances along the squeezed quadrature 
as Vs = 0.1, the entangling beam splitter is symmetric 
(Ts = 50%) and the tap beam splitter has a transmission 
of T = 0.7. The distillation with a pure entangled state 
(Va = 1/Ks) as well as a mixed state (V A = I/V5 + 10) 
are investigated as a function of the success probability, 
and shown in Fig. [3] (a) and (b) respectively. Follow- 
ing the theory of section II A 1[ we calculate the upper 
bound on distillable entanglement of the non-Gaussian 
state before distillation as shown by the bold straight 
lines in Fig. [3j The lower bounds on distillable entangle- 
ment after distillation are computed and shown in Fig. [3] 
by the dashed lines. We see that the proof of entangle- 
ment distillation of non- Gaussian states already fails for a 
mixed state with a small amount of excess noise. The up- 
per bounds on Gaussian entanglement after distillation 
are also computed and shown in Fig. [3] by solid lines. 
As the success probability reduces, the Gaussian entan- 
glement increases. Furthermore, when it surpasses the 
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upper bound on distillable entanglement before distilla- 
tion (bold solid lines) at a certain low success probability, 
the distilled state is Gaussified as well and thus we can 
justify a Gaussian state in the entanglement measure. 

(a) distillation with initially pure state 



h 2.5 




Success probability P s (log scale) 

(b) distillation with initially mixed states 
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Success probability P s (log scale) 

FIG. 3. (color online) Theoretical simulations of distillable en- 
tanglement of non-Gaussian mixed states as a function of suc- 
cess probability. The two plots are corresponding to two dif- 
ferent purities of the entangled input states. In figure (a), the 
distillation with initially pure state is plotted (Va = 1/Vs), 
while figure (b) shows the distillation with initially mixed 
states (Va — 1/Vs + 10). The other parameters are taken 
as: V s = 0.1, V N = 1, T s = 0.5, T = 0.7, m = 1, m = 0, 
pi = 0.2, p2 = 1 — pi. In both plots, the lower dashed line 
shows the lower bound on distillable entanglement, and the 
upper solid line is the upper bound on Gaussian entanglement. 
The bold straight line is the upper bound of the non- Gaussian 
distillable entanglement before distillation. 



III. EXPERIMENTAL REALIZATION 

The experimental realization of the distillation of 
corrupted entangled states consists of three parts as 



schematically illustrated in Fig. [4j the preparation, dis- 
tillation and verification. In the following we describe 
each part. 

A. Generation of polarisation squeezing and 
entanglement 

The generation of polarization squeezed beams serves 
as the first step for the demonstration of entanglement 
distillation. Here we exploit the Kerr nonlinearity of sil- 
ica fibers experienced by ultrashort laser pulses for the 
generation of quadrature squeezed states. Fig. [5] depicts 
the setup for the generation of a polarization squeezed 
beam. A pulsed (140 fs) Cr 4+ :YAG laser at a wave- 
length of 1500 nm and a repetition rate of 163 MHz is 
used to pump a polarization-maintaining fiber. Two lin- 
early polarized light pulses with identical intensities are 
traveling in single pass along the orthogonal polarization 
axes (x and y) of the fiber. Two quadrature squeezed 
states, the squeezed quadrature of which are skewed by 
# sq from the amplitude direction, are thereby indepen- 
dently generated. After the fiber the emerging pulses 
are overlapped with a tt/2 relative phase difference. The 
relative phase difference is achieved using a birefringence 
pre-compensation, an unbalanced Michelson-like interfer- 
ometer [38-41 j. This is controlled by a feedback lock- 
ing loop based on a S2 measurement of a small portion 
(<0.1%) of the fiber output. The measured error sig- 
nal is fed back to the piezo-electric element of the pre- 
compensation via a PI controller, so that the S2 parame- 
ter of the output mode vanishes. This results in a circu- 
larly polarized beam at the fiber output ((Si) = (S2) = 0, 
(S3) = (So) = a 2 ). The corresponding Stokes operator 
uncertainty relations are reduced to a single nontrivial 
one in the so-called S1—S2 dark plane: A 2 Sq A 2 Sq+ 7V /2 > 
\(S 3 )\ 2 , where S(0) = cos(#)Si + sin(0)S 2 denotes a gen- 
eral Stokes parameter rotated by in the dark Si — S2 
plane with (Sq) = 0. Therefore, polarization squeezing 
occurs if A 2 So < \(Ss)\ = a 2 , in which A 2 Sq can be di- 
rectly measured in a Stokes measurement [39 . As the 
noise of Stokes parameters Sq is linked to the quadra- 
ture noise of the Kerr squeezed modes in the same an- 
gle (A 2 S e « a(SX Xi9 - SX yt0 )/y/2 w a 2 A 2 X e [39 ), the 
squeezed Stokes operator is S(9 sq ) and the orthogonal, 
anti-squeezed Stokes operator is S(0 S(l + tt/2). Due to 
the equivalence between the polarization squeezing and 
vacuum squeezing [42], we utilize the conjugate quadra- 
tures X and P to denote the polarization squeezed and 
anti-squeezed Stokes operators. 

To generate polarization entanglement two identical 
polarization-maintaining fibers are used. Two polariza- 
tion squeezed beams, labeled A and B, are then gen- 
erated. By balancing the transmitted optical power of 
the two fibers, the two resultant polarization squeezed 
beams have identical squeezing angles, squeezing and 
anti-squeezing properties. The two polarization squeezed 
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FIG. 4. Schematics of the experimental setup for the preparation, distillation and verification of the distillation of entanglement 
from a non-Gaussian mixture of polarization entangled states. 
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FIG. 5. Setup for the generation of polarisation squeezing. 
The fiber is a 13.2 meters long polarization-maintaining 3M 
FS-PM-7811 fiber with a mode field diameter of 5.7 /im and 
a beat length of 1.67 mm. The interferometer in front of the 
fiber introduces a phase shift 5<j) between the two orthogonally 
polarised pulses to pre-compensate for the birefringence. A/4, 
A/ 2: quarter-, half-wave plates, PBS: polarising beam split- 
ter. PZT: Piezo-electric element. 



beams are then interfered on a 50/50 beam splitter 
(Fig. [3 with the interference visibility aligned to be 
> 98 The relative phase between the two input beams 
is locked to tt/2 so that the two output beams after the 
beam splitter have equal intensity and are maximally en- 
tangled. The two entangled outputs remain circularly 
polarised, thus the quantum correlations between them 
are lying in the dark Si — §2 plane with the signatures 
S A (0 sq ) + S B (0 sq ) and S A (6 sq + tt/2) - S B (6 sq + 
tt/2) -> (or X A + X B -> and P A - P B -> 0). 



0.9/44. The entangled beam is then transmitted through 
the lossy channel with 45 realizations. Combining all 
these realizations a non-Gaussian mixed state, such as 
the one described by eqn. (12), is achieved, with the 
probabilities pi all being identical. However, after the 
measurement we can select a certain probability enve- 
lope function to give the different channels pre-specified 
probability weights. With this technique we can eas- 
ily implement different transmission scenarios (see e.g. 
Fig. [9}l, Fig. [10} 1, and Fig. [ll}l). As a result of the 
lossy transmission, the Gaussian entanglement between 
the two beams A and B are degraded or completely lost. 



C. Entanglement distillation 

The distillation operation consists of a measurement of 
X on a small portion of the mixed entangled beam. This 
is implemented by tapping 7% of beam B after the ND 
filter using a beam splitter. The measurement is followed 
by a probabilistic heralding process where the remaining 
state is kept or discarded, conditioned on the measure- 
ment outcomes: e.g. if the outcome of the weak mea- 
surement is larger than the threshold value, Xth, then the 
state is kept. Note that the signal heralding process could 
in principle be implemented electro-optically to gener- 
ate a freely propagating distilled signal state. However, 
to avoid such complications, our conditioning is instead 
based on digital data post-selection using a verification 
measurement on the conjugate quadratures X and P of 
the beams A and B. 



B. Preparation of a non-Gaussian mixed state 

The preparation of a non- Gaussian mixture of polar- 
ization entangled states is implemented by transmitting 
one of the entangled beams, e.g. beam B, through a 
controllable neutral density filter (ND). The filter is im- 
plemented to produce a lossy channel with N = 45 differ- 
ent transmittance levels, ranging from 0.1 to 1 in steps of 



D. The tap and verification measurement 

The tap and verification measurement are accom- 
plished simultaneously by three independent Stokes mea- 
surement apparatuses. Each measurement apparatus 
consists of a half- wave plate and a polarizing beam split- 
ter (PBS). Since the light beam is circularly (S3) po- 
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larized, a rotation of the half- wave plate enables the 
measurement of different Stokes parameters lying in the 
'dark' polarization plane. For the tap measurement the 
half-wave plate is always set at the angle corresponding 
to X in the 'dark' polarization plane. Via the verifi- 
cation measurement setup, the Gaussian properties of 
the entangled states are characterized by measuring the 
entries of the covariance matrix. By generating near 
symmetric states and choosing a proper reference frame, 
we assume that the intra-correlations (such as (XaPa)) 
are zero. The measurements of these entries are accom- 
plished by applying polarization measurements of beam 
A and B with both the half-wave plates set to the an- 
gle corresponding to either X or P in the 'dark' polar- 
ization plane. The outputs of the PBS are detected by 
identical pairs of balanced photo-detectors based on 98% 
quantum efficiency InGaAs PIN-photodiodes and with an 
incorporated low-pass filter in order to avoid ac satura- 
tion due to the laser repetition oscillation. The detected 
AC photocurrents are passively pairwise subtracted and 
subsequently down-mixed at 17 MHz, low-pass filtered 
(1.9 MHz), and amplified (FEMTO DHPVA-100) before 
being oversampled by a 16-bit A/D card (Gage CompuS- 
cope 1610) at the rate of 10 7 samples per second. The 
time series data are then low-passed with a digital top- 
hat filter with a bandwidth of 1 MHz. After these data 
processing steps, the noise statistics of the Stokes param- 
eters are characterized at 17 MHz relative to the optical 
field carrier frequency (^200 THz) with a bandwidth of 
1 MHz. The signal is sampled around this sideband to 
avoid classical noise present in the frequency band around 
the carrier [44]. For each polarization measurement, the 
detected photocurrent noise of beam A and B and the tap 
beam were simultaneously sampled for 2.4 x 10 8 times, 
thus the self and cross correlations between the data set 
of A and B could thereby easily be characterized. The 
covariance matrix was subsequently determined and the 
log-negativity was calculated according to eqn. (|7|). 

IV. EXPERIMENTAL RESULTS 

For perfect transmission (corresponding to no loss in 
beam B), the marginal distributions of the entangled 
beams, A and B, along the quadrature X and P are plot- 
ted in Fig. [6j In Fig. [6ja) the procedure of realizing the 
different noisy channels is shown. The sampled data of 
different attenuation channels is concatenated according 
to the different weights of the transmission probabilities. 
These samples then provide the measurement data for 
the distillation procedure. From Fig. [6^a) and Fig. [6] (b) 
we can see that, each individual mode exhibits a large 
excess noise (measured fluctuation > 17 dB). However 
the joint measurements on the entangled beams A and 
B exhibit less noise fluctuation than the shot noise ref- 
erence, as shown in Fig. [6^c). The observed two- mode 
squeezing between beam A and B is — 2.6 ± 0.3 dB and 
—2.4 ±0.3 dB for X and P, respectively. From the deter- 



mined covariance matrix we compute the log-negativity 
to be 0.76 ± 0.08. 

a) Measurement on beam A 
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FIG. 6. (color online). Experimentally measured marginal 
distributions associated with the (a) X and P of beam A, (b) 
X and P of beam B and (c) the joint measurements Xa + Xb 
and Pa — Pb • The black and red curves are the distributions 
for shot noise and the quadrature on measurement, respec- 
tively. 

To experimentally demonstrate the distillation of en- 
tanglement out of non-Gaussian noise, three different 
lossy channels are considered: the discrete erasure chan- 
nel, where the transmission randomly alternates between 
two different levels, and two semi-continuous channels, 
where the transmittance alternates between 45 different 
levels with specified probability amplitudes. The proba- 
bility distributions of the transmittance for the discrete 
channel and the continuous channels are shown in Fig. [9} 
1, Fig. [TO]- 1, and Fig. [TTJ-1, respectively. 



A. The discrete lossy channel 

The discrete erasure channel alternates between full 
(100%) transmission and 25% transmission at a proba- 
bility of 0.5. Each realisation is concatenated to each 
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other with identical weights. The concatenation proce- 
dure yields the same statistical values as true randomly 
varying data. After transmission the resulting state is 
a mixture of a highly entangled state and a weakly en- 
tangled state. In the inset of Fig. [7J we show marginal 
distributions illustrating the single beam statistics of the 
individual components of the mixture. The statistics of 
beam B is seen to be contaminated with the attenuated 
entangled state thus producing non- Gaussian statistics. 
For this state we measure the correlations in X and P to 
be above the shot noise level by 5.5 ±0.3 dB and 5.6 ±0.3 
dB, respectively, and the Gaussian LN to be — 1.63±0.02. 
The Gaussian entanglement is completely lost as a result 
of the introduction of such time-dependent loss. This 
is in stark contrast to the scenario where only station- 
ary loss (corresponding to Gaussian loss) is inflicting the 
entangled states. In that case, a certain degree of Gaus- 
sian entanglement will always survive, although it will be 
small for high loss levels. 

The state is then fed into the distiller and we perform 
homodyne measurements on beam A, beam B and the 
tap beam simultaneously. By measuring X in the tap we 
construct the distribution shown by the red curve in the 
left hand side of Fig. (7^a). The data trace of the mixed 
tap signal is plotted accordingly on the right hand side. 
The measurements of X and P of the signal entangled 
states were recorded as well. For simplicity, we only show 
the distribution for the beam B (in Fig.[7Jb-l),(b-2)) and 
the joint distribution of beam A and B (in Fig. [7jc-l),(c- 
2)). The blue (dashed) and red curves denote the dis- 
tributions before and after the post-selection process, re- 
spectively. From the blue curves shown in Fig. [7| we can 
see that the entanglement between A and B is lost due 
to the non-Gaussian noise. Performing postselection on 
this data by conditioning it on the tap measurement out- 
come (denoted by X t h = 9-0), we observe a recovery of 
the entanglement. That is, the correlated distribution of 
the signal turns out to be narrower than that of the shot 
noise (as shown by red curves in Fig. [TJc). 

Using the data shown in Fig. [7J a tomographic recon- 
struction of the covariance matrices of the distilled en- 
tangled state was carried out. From these data we de- 
termined the most significant eight of the ten indepen- 
dent parameters of the covariance matrix, namely the 
variances of four quadratures P^, Pb and co- 

variances between all pairs of quadratures of the entan- 
gled beams A and B. As mentioned before, the intra- 
correlations were ignored. The resulting covariance ma- 
trices are plotted in Fig. [8] for ten different postseletion 
threshold values from X t h = 0.0 to X t h = 9.0 with a step 
of 1.0. With increasing postselection threshold the distil- 
lation becomes stronger, as shown by the reduction of the 
quadrature variances of Xa-> Xb and the increase of the 
quadrature variances of Pa, Pb- Moreover, the reduction 
(or increase) of the covariances C(Xa^b) (C(Pa, Pb)) 
was shown slightly slower. Consequently, the entangle- 
ment of the two modes A and B was enhanced by the 
distillation. 
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FIG. 7. (color online). Experimentally measured marginal 
distributions illustrating the effect of distillation, (a) Exam- 
ple of concatenated sampled data and the resulting marginal 
distribution for the amplitude quadrature in the tap measure- 
ment. The vertical line indicates the threshold value chosen 
for this realization, (b) Marginal distributions associated with 
the measurements of X and P of beam B (two left figures) 
and (c) the joint measurements Xa + Xb and Pa — Pb (two 
right figures). The black, blue and red curves are the distri- 
butions for shot noise, the mixed state before distillation and 
after distillation, respectively. Inset: phase-space representa- 
tion of the non-Gaussian mixed state and the post-selection 
procedure used in the measurements. The black vertical line 
indicates the threshold value. 

Furthermore, the distilled entanglement, or log- 
negativity, was investigated as a function of the suc- 
cess probability, as shown in Fig. [9] by black open cir- 
cles. The error bars of the distilled log- negativity de- 
pend on two contributions: First the measurement er- 
ror, which is mainly associated with the finite resolu- 
tion of the A/D converter and noise of the electronic 
amplifiers. This is considered by estimating the exper- 
imental error for all the elements of the covariance ma- 
trix as '0.03'. The measurement error for the LN can 
be simulated by a Monte-Carlo model. Second, the sta- 
tistical error is due to the finite measurement time and 
the postselecti on process . It is considered by adding a 
scaled term ^/2/(N — 1), where N denotes the number 
of postselected data [45]. The probability distributions 



11 




FIG. 8. (color online). Reconstructed covariance matrices of 
distilled entangled states. The brown segmented plane shows 
the region for the individual elements in the covariance ma- 
trix. The sub-bars represent the results of our distillation 
protocol for 10 different threshold values postseletion thresh- 
old values from Xth = 0.0 to Xth = 9.0 with a step of 1.0. 



of the two superimposed states in the mixture after dis- 
tillation are shown for different post select ion thresholds, 
corresponding to X t h = 0.0,2.0,4.0,6.0,9.0, labeled by 
1-5 in order. The plots explicitly show the effect of the 
distillation protocol, when the postselection threshold in- 
creases, the Gaussian LN increases, ultimately approach- 
ing the LN of the input entanglement without losses. The 
probability distribution tends to a single valued distribu- 
tion, therefore the mixture of the two Gaussian entan- 
gled states reduces to a single highly entangled Gaussian 
state, thus demonstrating the act of Gaussification. How- 
ever, the amount of distilled data, or success probability, 
decreases, causing an increase in the statistical error on 
the distillable entanglement. Based on the experimental 
parameters, a theoretical simulation is plotted by the red 
curve and shows a very good agreement with the exper- 
imental results. 

To further investigate whether the total entanglement 
is increased after distillation, we compute the upper 
bound for the LN before distillation and verify that this 
bound can be surpassed by the Gaussian LN after distil- 
lation. The upper bound of LN without the Gaussian ap- 
proximation is computable from the LN of each Gaussian 
state in the mixture [34 , and we find LN upper = 0.49, 
which is shown in Fig.[9]by the dashed black line. We see 
that for a success probability around 10 -4 the Gaussian 
LN crosses the upper bound for entanglement. Since the 
state at this point is perfectly Gaussified we may con- 
clude that the total entanglement of the state has in- 
deed increased as a result of the distillation. Fig. |9}5 
gives another explicit explanation by showing that the 
probability contribution from the 75% attenuated data 
reaches when the post selection threshold is set to 
X t h = 9.0, which corresponds to the distilled entangle- 
ment of LNg = 0.67 ± 0.08 with a success probability 
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FIG. 9. (color online). Experimental and theoretical results 
outlining the distillation of an entangled state from a discrete 
lossy channel. The experimental results are marked by cir- 
cles and the theoretical prediction is plotted by the red solid 
line. The bound for Gaussian entanglement is given by the 
blue line, and the upper bound for total entanglement before 
distillation is given by the black dashed line. Both bounds 
are surpassed by the experimental data. The weight of the 
two constituents in the mixed state after distillation for var- 
ious threshold values is also experimentally investigated and 
shown in the plots labeled by 1-5. The error bars of the log- 
negativity represent the standard deviations. 



P s = 1.69 x 10~ 5 . On the other hand, from Fig. [^3 and 
Fig. [9]4, we see that even a small contribution from the 
75% attenuated data will reduce the useful entanglement 
for Gaussian operations. 



B. The continuous lossy channel 

We now generalize the lossy channels to have a con- 
tinuously transmittance distribution. The channel trans- 
mittance distribution is simulated by taking 45 differ- 
ent transmission levels as opposed to the two levels in 
the previous section. In Ref. [27 we reported a chan- 
nel whose transmittance is given by an exponentially de- 
caying function with a long tail of low transmittances, 
which simulates a short-term free-space optical commu- 
nications channel where atmospheric turbulence causes 
scattering and beam pointing noise [46] . We showed that 
the entanglement available for Gaussian operations can 
be successfully distilled from -0.11 ± 0.05 to 0.39 ± 0.07 
with a success probability of 1.66 x 10 -5 . However, in 
practical scenarios for a transmission channel, the highest 
transmittance level may not have the biggest weight in 
the probability distribution and therefore the distributed 
peak may be displaced from the 100% transmittance 
level. Further, there might be more than one peak in 
the probability distribution diagram. For instance, due 
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to some strong beam pointing noise another distributed 
peak will appear in the area of low transmittance lev- 
els. In the following we will test the performance of the 
distillation protocol for two different transmittance dis- 
tributions. First, when the mixed state has a peak of the 
transmittance distribution which is displaced from 1 to 
0.8 (Fig. [To]-!). Second, when we incorporate a second 
peak which is located around the transmittance level of 
0.3 (Fig.[ll}l). 
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FIG. 10. (color online). Experimental and theoretical re- 
sults outlining the distillation of an entangled state from a 
simulated continuous lossy channel in which the peak of the 
transmittance distribution is displaced from 1 to 0.8. The 
experimental results are marked by circles and the theoret- 
ical prediction by the red solid curve. The evolution of the 
weights of the various constituents in the mixed state as the 
threshold value is changed is shown in the figures labeled 1-5. 
The error bars of the log- negativity represent the standard 
deviations. 

As shown in Fig. [TUJ after propagation through the 
one-peak displaced channel the Gaussian LN of the mixed 
state is found to be —0.50 ± 0.04, which is below the 
bound for available entanglement (shown by the solid blue 
line) and substantially lower than the original value of 
0.76 ± 0.08. The state is subsequently distilled and the 
change in the Gaussian LN as the threshold value in- 
creases (and the success probability decreases) was in- 
vestigated both experimentally (black open circles) and 
theoretically (red curve). The evolution of the mixture 
is directly visualized in the series of probability distribu- 
tions in Fig. [10]T to[l0}5 corresponding to the postselec- 
tion thresholds X^=0.0, 3.0, 5.0, 7.0, 9.0 respectively. 
We see that the distribution weights of the low trans- 
mittance levels is gradually reduced, while the weights 
of the high transmittance levels is increased as the post- 
selection process becomes more and more restrictive by 
increasing the threshold value. E.g. for X t h = 9 the 
probabilities associated with transmission levels lower 
than 0.7 are decreased from 20% before distillation to 



1.4% and the probability for transmission levels higher 
than 0.7 transmission are increased to 98.6% as opposed 
to 80% before distillation. It is thus clear from these 
figures that the highly entangled states in the mixture 
have larger weight after distillation, and the correspond- 
ing Gaussian LN after distillation rises to 0.19 ±0.06 with 
the success probability of 5.16 x 10 -5 . 
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FIG. 11. (color online). Experimental and theoretical results 
outlining the distillation of an entangled state from a sim- 
ulated continuous lossy channel in which the transmittance 
levels are distributed as such that there are two peaks at both 
high transmittance levels (0.8) and low transmittance levels 
(0.3). The experimental results are marked by circles and the 
theoretical prediction by the red solid curve. The evolution of 
the weights of the various constituents in the mixed state as 
the threshold value is changed is shown in the figures labeled 
1-5. The error bars of the log-negativity represent standard 
deviations. 

We now turn to investigate the distillation after propa- 
gation through the two-peak displaced channel as shown 
in Fig. [TTJ-1. Before distillation the Gaussian LN of 
the mixed state is found to be — 1.13 ± 0.02. Like- 
wise, the relation between the distilled Gaussian LN and 
the success probability was investigated both experimen- 
tally and theoretically. The results are shown in Fig. [TT] 
by black open circles and the red curve, respectively. 
Through the probability distribution plots in Fig. [TT] -1 
to[TT]-5, the evolution of the mixture corresponding to dif- 
ferent choices of postselection thresholds (X^=0.0, 3.0, 
5.0, 7.0, 9.0 respectively) was illustrated with the same 
trend that we see on the distillation after the one-peak 
displaced channel. For X t h = 9 the probabilities associ- 
ated with transmission levels lower than 0.7 are decreased 
from 48% before distillation to 1.6% and the probability 
for transmission levels higher than 0.7 transmission are 
increased to 98.4% as opposed to 52% before distilla- 
tion, and the corresponding Gaussian LN after distilla- 
tion reaches 0.19 ± 0.06 with the success probability of 
3.39 x 10- 5 . 
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After having shown the successful entanglement dis- 
tillation on different distributions of non-Gaussian noise, 
we should note that the successful entanglement distil- 
lation depends on the transmittance distribution of the 
lossy channel. For some distributions, the success prob- 
ability for distilling available entanglement for Gaussian 
operations will be extremely small or not be possible. 
For example, after a channel with the transmittance uni- 
formally distributed, the Gaussian log-negativity LN$ = 
— 1.26 ± 0.02 before distillation will only be increased to 
-0.76±0.03 with a success probability of 1.32 x 10" 5 . In 
general more uniform transmittance distributions turned 
out to be more difficult for the distillation procedure. 
Distributions with high probabilities for high transmis- 
sion levels and pronounced tails and peaks at low trans- 
mission levels (as would be expected in atmospheric chan- 
nels) are more suited. 

V. SUMMARY 

In summary, we have proposed a simple method of dis- 
tilling entanglement from single copies of quantum states 



that have undergone attenuation in a lossy channel with 
varying transmission. Simply by implementing a weak 
measurement based on a beam splitter and a homodyne 
detector, it is possible to distill a set of highly entan- 
gled states from a larger set of unentangled states if the 
mixed state is non-Gaussian. The protocol was success- 
fully demonstrated for a discrete erasure channel where 
the transmittance alternates between 2 levels and two 
semi-continuous transmission channels where the trans- 
mission levels span 45 levels with specified distributions, 
respectively. We show that the degree of Gaussian en- 
tanglement (which is relevant for Gaussian information 
processing) was substantially increased by the action of 
distillation. Moreover, we proved experimentally that the 
total entanglement was indeed increased for the discrete 
channel. We found that the successful entanglement dis- 
tillation depends on the transmittance distribution of the 
lossy channel. The demonstration of a distillation proto- 
col for non-Gaussian noise provides a crucial step towards 
the construction of a quantum repeater for transmitting 
continuous variables quantum states over long distances 
in channels inflicted by non-Gaussian noise. 
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